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Abstract.  A  new  computational  method  for  characterizing  the  relationship  between  surface  crystallography  and 
electrical  conductivity  in  anisotropic  materials  with  lowT  frequency  electromagnetic  techniques  is  presented.  The 
method  is  discussed  from  the  standpoint  of  characterizing  the  orientation  of  a  single  grain,  as  well  as  characterizing 
statistical  information  about  gram  ensembles  in  the  nncrostructure.  Large-area  election  backscatter  diffraction  (EBSD) 
data  was  obtained  and  used  in  conjunction  with  a  synthetic  aperture  approach  to  simulate  the  eddy  current  response  of 
beta  annealed  Ti-6A1-4V.  Experimental  eddy  cuirent  results  are  compared  to  the  computed  eddy  current 
approximations  based  on  electron  backscatter  diffraction  (EBSD)  data,  demonstrating  good  agreement  The 
detectability  of  notches  in  the  presence  of  noise  from  microstructure  is  analyzed  with  the  described  simulation  method 
and  advantages  and  limitations  of  this  method  are  discussed  relative  to  other  NDE  techniques  for  such  analysis. 

Keywords:  Eddy  Cuirent.  Microstructure.  Hexagonal.  Anisotropic 
PACS:  81.70.-q.  81.05.Xj 


INTRODUCTION 

One  of  the  main  confounding  factors  in  using  the  methods  of  eddy  current  testing  (ECT)  to  detect  small  (<  0.1 
mm)  cracks  in  engineering  anisotropic  alloys  is  the  microstructure  of  the  material.  When  the  change  of  impedance 
due  to  the  presence  of  the  crack  is  on  the  same  order  as  the  change  due  to  perturbations  in  crystallographic  orientation, 
the  crack  signal  gets  lost  in  what  is  often  perceived  as  the  surrounding  "grain  noise”  [1],  This  situation  is  illustrated 
in  Fig  1.  The  same  sized  electrical  discharge  machined  (EDM)  notch  was  scanned  in  two  different  materials.  The 
main  difference  between  the  two  images  is  that  the  material  in  the  image  on  the  left  is  isotropic  vs.  the  anisotropic 
material  on  the  right.  Several  methods  could  be  used  to  enhance  the  signal  from  the  notch  in  Fig.  1(b).  Any  of  several 
image  processing  routines  could  be  used  to  identify  and  subtract  the  features  caused  by  the  microstructure  from  that 
due  to  the  presence  of  damage.  This  would  offer  some  sensitivity  to  the  notch  or  it  could  provide  a  method  to  reduce 
the  noise.  This  would  be  a  similar  approach  to  using  spatial  Fourier  transforms  to  extract  features  from  miages.  These 
methods  do  offer  some  enhanced  sensitivity  to  the  notch  signal,  but  in  the  process  they  make  parameter  estimation 
much  more  involved.  For  NDE  purposes,  the  ideal  situation  is  one  in  which  this  data  can  be  used  in  its  raw  form  to 
estimate  characteristics  of  the  damage.  For  this  type  of  analysis,  a  thorough  understanding  of  the  physics  involved  is 
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FIGURE  1.  Eddy  current  scans  taken  over  the  same  sized  crack  in  (a)  A L  2024  and  in  (b)  Ti-6A1-4V  (P  annealed)  [1]. 
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FIGURE  2.  linage  of  the  sample  analyzed  in  this  work  The  grains  are  clearly  large  enough  that  the  coil  can  resolve  thein. 

extremely  helpful.  To  this  end.  this  work  is  mainly  concerned  with  proof  of  the  theory  that  this  noise  is  related  to 
micro  structure  and  in  providing  a  potential  forward  model  approximation  for  this  type  of  problem. 

EXPERIMENTS 


Experiments  were  performed  initially  to  determine  the  source  of  the  material  noise  and  then  to  provide  sample 
information  to  the  forward  models  for  verification.  An  appropriate  sample  was  chosen  for  the  purposes  of  the  work, 
and  then  several  different  experimental  techniques  were  used  to  collect  information  about  the  microstructure.  The 
different  experimental  data  sets  were  compared  for  verification. 

Sample 

In  order  to  illustrate  and  exacerbate  the  role  of  gram  orientation  on  conductivity,  thereby  demonstrating  the  source 
of  microstructure-induced  noise,  a  coarse-grained  beta-annealed  Ti-6A1-4V  sample  was  selected  for  this  work.  The 
grains  were  sufficiently  large  that  the  diameter  of  a  standard,  commercially  available  ECT  probe  was  less  than  the 
typical  grain  diameter  and  therefore  the  signal  at  a  given  could  be  due  to  a  single  grain.  It  is  noteworthy  that  the  large 
grain  dimensions,  up  to  -  1  cm.  in  this  sample  imply  a  large  radius  of  curvature  of  the  colony  boundaries,  much  larger, 
in  fact,  than  the  depth  contributing  to  the  eddy  current  signal  and  thus  complications  due  to  subsurface  variations  in 
orientation  that  may  be  present  in  fine  grained  materials  were  not  a  source  of  variability  here.  The  optical  images  of 
the  sample  shown  in  Fig.  2  illustrate  this  point  more  clearly.  T i-6 A1-4V  consists  primarily  of  the  hexagonal  close 
packed  a-phase  (-93%  volume  fraction)  with  the  remainder  balanced  by  the  body  center  cubic  (3-phase.  The  a-phase 
exhibits  isotropy  in  the  basal  plane,  but  most  physical  properties  are  moderately  anisotropic  and  depend  on  angle 
relative  to  the  [0001]  (c-axis)  direction.  The  variation  in  conductivity  from  the  basal  plane  to  the  c-axis  is  nominally 
about  6%  difference,  and  the  conductivity'  in  the  basal  plane  is  assumed  to  be  on  the  order  of  2.4rlOG  S/m  [2], 

Electron  Backscatter  Diffraction 


Since  most  properties  vary  as  a  function  of  orientation,  and  the  cost  of  these  measurements  is  becoming  more 
affordable,  large-area  microstructure  characterization  involving  measurement  of  local  orientation  via  EBSD  is 
becoming  increasingly  commonplace  (e  g.  [3]).  Crystallographic  orientation  in  an  orthonormal  basis  can  be  fully 


FIGURE  3.  (a)  The  hexagonal  crystal  in  an  orientation  in  line  with  the  sample  axes,  (b)  rotated  an  angle  <pi  about  the  initial 
z-axis.  (c)  rotated  an  angle  6  about  the  new  x-axis.  and  (d)  rotated  an  angle  y>;  about  the  newest  z-axis. 


2 

Distribution  A.  Approved  for  public  release;  distribution  unlimited. 


FIGURE  4. 


Crystal  orientation  map  of  the  coarse-grained  Ti-6A1-4V  sample.  Here,  color  corresponds  to  crystallographic 
orientation.  The  large  grain  in  the  center  is  outlined  for  reference. 


described  by  a  set  of  Euler  angles  ((pv  6,  cpz)  in  Bunge's  active  form,  as  shown  in  Fig  3.  The  Euler  angles  essentially 
define  the  rotations  necessary'  to  bring  the  sample  axes  into  coincidence  with  an  arbitrarily  oriented  crystal  [4].  In 
this  work,  a  direct  comparison  between  eddy  current  images  and  the  actual  sample  cannot  be  made  without  this 
orientation  information.  Orientation  data  were  collected  with  electron  backscatter  diffraction  using  a  Phillips  XL30 
scanning  electron  microscope  (SEM)  operated  at  an  accelerating  voltage  of  20  kV  at  a  probe  current  of  50  nA.  The 
resulting  data  is  shown  in  Fig  4. 


Eddy  Current  Experiments 

A  commercially  available.  2MHz  center  frequency  probe  was  used  to  collect  EC'T  data  from  the  large  grained 
sample  shown  above.  The  probe  was  connected  to  a  Nortec  2000D  eddy  current  scope  and  raster  scanned  over  the 
entire  specimen.  This  combination  of  instrument  and  probe  has  been  shown  to  produce  detectable  signals  from 
conductivity  changes  of  below  6%.  which  is  the  maximum  change  that  would  be  expected  from  the  sample  (i.e.  the 
percent  difference  between  the  c-axis  and  the  basal  plane  conductivities).  To  achieve  such  sensitive  measurements, 
one  crucial  step  is  to  determine  the  probe’s  sensitivity  to  conductivity  vs.  liftoff.  A  calibration  curve  for  this  probe 
generated  with  physical  models  [5]  can  be  seen  in  Fig.  5.  The  conductivity  range  that  the  sample  is  in  is  close  to  1 .75 
%IACS,  which,  upon  inspection  of  the  curve  in  that  range,  means  that  liftoff  is  not  orthogonal  to  conductivity  changes. 
However,  with  conductivity  changes  on  the  order  of  6%  [2],  the  magnitude  change  is  close  enough  to  linear  to  use  it 
as  the  impedance  measurement.  Liftoff  was  controlled  by  using  the  probe  in  contact  with  the  specimen.  An  image  of 
the  magnitude  of  impedance  above  the  sample  is  shown  m  Fig.  6.  This  image  shows  some  of  the  grains  that  can  be 
seen  in  the  OIM  image  indicating  that  the  ECT  method  is  capable  of  distinguishing  the  grains,  and  that  the  change  in 
the  impedance  of  the  coil  is  due  mainly  to  the  conductivity  change  from  one  region  to  another. 

APPROXIMATION  METHOD 

Knowing  that  this  variation  in  impedance  measurements  over  Ti  and  its  alloys  is  a  conductivity  variation  due  to 
microstructure,  a  quick  and  accurate  approximation  method  for  forward  models  of  this  phenomenon  can  be  obtained. 
This  is  useful  if  statistical  quantification  of  the  microstrucnire  is  to  be  accomplished  with  ECT  methods.  The  key  to 
the  approximation  lies  in  several  assumptions  that  will  be  enumerated  here  to  clarify-  the  method  as  well  as  provide 
potential  avenues  to  refine  the  approximation. 


r-pedarce  Pane  fw  yx  2MHz  Cd 


FIGURE  5.  Calibration  curve  for  the  eddy  cunent  probe.  The  curve  represents  change  in  conductivity,  and  the  hash  marks 

represents  a  change  m  liftoff  at  a  fixed  conductivity. 
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FIGURE  6.  ECT  scan  of  the  sample. 

Assumption  1:  Uniform  Electric  Field 

This  assumption  states  that  underneath  the  coil,  the  electric  field  is  uniform  inside  the  conductor.  This  is  an 
assumption  since,  as  can  be  seen  in  any  mathematical  analysis  of  eddy  current  problems  [6.  7],  the  electric  field  m  the 
conductor  is  dependent  on  the  conductivity  of  the  sample.  However,  this  assumption  hinges  on  the  fact  that  the 
conductivity  changes  are  so  small  that  the  variation  of  electric  field  from  one  grain  to  the  other  is  small.  The  definition 
of ‘  small  enough"  here  is  not  analyzed,  and  in  fact  this  assumption  might  be  fairly  stringent.  However,  since  the  goal 
is  to  produce  an  approximation  of  an  eddy  current  scan,  this  assumption  will  be  enforced. 

Assumption  2:  No  Changes  in  Depth 

The  second  assumption  is  that  the  grains  that  we  see  in  the  EBSD  map  are  ‘pillars",  or  that  they  are  constant  in  the 
depth  of  the  sample  This  is  clearly  a  strict  and  rather  limiting  assumption.  The  skin  depth  of  the  eddy  currents  in  this 
material  at  2MHz  is  around  230pm  Common  commercially  used  Ti  alloys  have  microstructures  with  grain  sizes  in 
the  range  of  10-20|im  meaning  that  the  present  analysis  would  be  confounded  by  non  surface-connected  grains  which 
would  contribute  to  the  eddy  current  response.  In  tins  case,  however,  the  large  grain  size  and  commensurate  large 
grain  boundary  curvature  implies  that  the  region  of  constant  orientation  far  exceeds  the  depth  of  the  eddy  currents. 
Optical  inspection  of  both  sides  of  the  sample  indicated  that  many  of  the  grains  actually  spanned  the  entire  thickness 
(~5  mm)  validating  that  assumption  that  the  EBSD  measurements  are  representative  of  the  orientations  at  least. 

Details  of  the  Method 

With  the  assumptions  clear,  the  approximation  routine  can  now  be  elaborated  upon.  The  first  thing  that  must  be 
understood  is  that  the  conductivity  in  the  material  is  different  for  electric  fields  of  different  direction.  The  constitutive 
relation  for  induced  current  density  in  the  presence  of  an  electric  field  is  [8]: 

/  =  (1) 

/  is  the  current  density,  a  is  the  conductivity  tensor,  and  E  is  the  electric  field.  If  we  assume  an  electric  field  in  the 
direction  of  a  unit  direction  vector,  r,  and  assign  it  some  arbitrary  magnitude,  say  Emag .  then  with  some  steps  of  linear 
algebra,  the  magnitude  of  the  current  density  induced  can  be  expressed  as: 

Jmag  =  Emagir*  I")1  2  (2) 

This  expresses  the  magnitude  of  the  current  density  as  a  function  of  an  electric  field  in  an  arbitrary  direction  and  the 
conductivity  tensor  of  the  sample.  The  current  density  here  is  not  necessarily  in  the  direction  of  the  electric  field,  but 
its  magnitude  is  correct.  We  can  define  an  effective  conductivity  as  aeff  =  (rrcr2r)1/2.  When  the  electric  field  is  in 
the  direction  of  one  of  the  principal  axes  of  the  conductivity  tensor,  the  conductivity  comes  out  as  the  single  crystal 
conductivity  in  that  direction,  as  expected.  Now.  for  an  arbitrary  set  of  rotations  given  Euler  angles,  a  rotation  matrix 
for  the  conductivity  tensor  can  be  defined  as  [4] : 
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r 


(3) 


R  = 


clci  -  c2sls3 
c3sl+clc2s3 

S2S3 


-cts3  -  c2c3s3 
C2C2C3  -  s3s3 
C3S2 


sls2 

~cls2 

c2 


C  represents  cosine  and  s  represents  sine,  and  the  subscripts  1.  2,  and  3  represent  the  operations  on  Euler  angles  <pL, 
Q.  and  <p~  respectively.  Given  a  new  direction.  Rr.  for  the  electric  field,  the  electric  field  can  be  expressed  as: 

E  —  EmagRr  (4) 

Again,  we  write  the  constitutive  relation  for  the  newly  aligned  region  as: 


/=  oE  =  Emag<rRr  (5) 

We  can  carry  out  the  same  operations  for  an  electric  field  in  an  arbitrary  direction  and  arrive  at  a  final  expression  for 
the  magnitude  of  the  current  density'  as: 


Jmag  =  Emag{rT  RT  t^Rr)12  (6) 

where  we’ve  used  the  fact  that  R~l  =  RT  The  effective  conductivity  in  this  expression: 

a -tff  =  (rT  RTo2Rr)1  2  (7) 

represents  the  conductivity  in  the  direction  of  the  electric  field. 

With  this  definition  of  effective  conductivity,  an  expression  for  the  conductivity  of  a  point  below  the  coil  can  be 
derived.  The  first  assumption  implies  that  the  field  below'  the  coil  and  m  the  sample  flow's  tangential  to  the  coil.  This 
is  a  similar  assumption  to  that  made  when  solving  an  axisymmetric  problem  with  scalar  conductivities  as  in  [5],  With 
that  in  mind,  the  direction  of  the  electric  field  at  a  point,  p,.  below'  the  coil  can  be  calculated  by  drawing  a  line  to  the 
point  and  then  drawing  the  normal  line  to  that  point  in  the  direction  opposite  the  flow  of  currents  in  the  coil  (due  to 
Lentz's  lau\  though  the  direction  doesn't  matter  m  this  application).  The  schematic  of  this  is  shown  in  Fig.  7.  Setting 
the  coil  center  as  the  origin  (0.0)  and  calling  the  vector  to  the  point  du  the  equation  for  the  tangential  direction.  r„ 
can  be  written  as: 


[0  0  1]  x  d{  [~dy  dx  0] 

- Wl - = - Mi -  * 

This  value  of  rh  along  with  the  Euler  angles  at  the  point  are  then  use  in  equation  7  to  calculate  the  effective 
conductivity  at  that  point.  Since  the  OIM  data  is  a  regular  grid  of  points,  averaging  every'  point  below  the  coil  gives 
what  is  similar  to  an  area-weighted  average  of  orientation-dependent  conductivities.  This  average,  w'hen  plotted  for 
every  coil  position,  gives  an  estimate  of  what  an  ECT  scan  would  look  like  above  an  anisotropic  material.  This 
essentially  amounts  to  a  synthetic  aperture  approach  for  eddy  current  scans  above  a  conductive  medium  with  subtle 
variations  in  conductivity. 


FIGURE  7.  The  direction  of  the  electric  field  used  in  effective  conductivity  calculations  at  a  point,  pt  below'  the  coil. 
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FIGURE  8.  Images  of  (a)  an  actual  ECT  scan  of  the  sample  and  (b)  a  simulated  ECT  scan. 


COMPARISONS  AND  VERIFICATION 


To  test  the  synthetic  aperture  routine,  the  Euler  angles  that  were  collected  from  the  sample  were  fed  into  the 
numerical  method  and  a  simulated  eddy  current  scan  was  performed.  The  images  are  shown  back-to-back  in  Fig  8. 
There  are  some  features  that  are  clearly  visible  in  both  images,  but  other  regions  are  very  difficult  to  compare.  The 
cause  of  this  is  partly  the  fact  that  the  first  assumption  implies  that  the  resxilts  of  the  simulation  are  not  physically 
accurate,  and  that  the  3D  grain  morphology  is  missing  from  the  simulation.  Again,  since  the  grains  are  mostly  large 
enough,  the  second  assumption  is  not  as  crucial,  but  it  could  be  relevant  if  the  angle  of  the  grain  boundary  is  drastic 
in  depth.  Further  validation  to  estimate  the  effect  of  the  assumptions  would  require  3D  sectioning  experiments  and 
the  use  of  a  full  numerical  solution  such  as  finite  element  method  (FEM).  boundary  element  method  (BEM),  or  volume 
integral  method  (VIM).  As  of  writing  this  paper,  these  experiments  have  not  been  performed. 

On  the  other  hand,  this  type  of  verification  has  been  demonstrated  to  a  lesser  extent  with  FEM  methods.  A  problem 
with  the  geometry  and  conductivities  shown  in  Fig.  9  was  analyzed  with  COMSOL  Multiphysics.  The  details  of  this 
model  can  be  found  in  [9],  The  two  different  conductivity  regions  wrere  simulated  using  a  different  rotation  matrix 
for  each.  The  FEM  simulation  gives  changes  in  impedance  of  the  coil,  whereas  the  synthetic  aperture  routine  give 
changes  in  conductivity.  Howrever.  it  was  already  shown  that  this  relationship  can  be  considered  lmear  with  such 
small  changes  in  conductivity,  so  the  results  from  each  were  normalized  and  plotted  in  Fig.  10.  The  synthetic  aperture 
routine  shows  very  good  agreement  with  the  FEM  results. 

CONCLUSIONS,  DISCUSSION,  AND  FUTURE  WORK 

From  the  data  provided  in  this  paper,  it  is  reasonable  to  conclude  that  noise  in  the  scans  shown  is  likely  due  to 
microstnicture  changes  in  an  anisotropic  material.  This  conclusion  has  been  reach  by  other  authors,  but  no  direct 
comparison  between  OIM  and  eddy  current  scans  has  been  made.  The  OIM  maps  of  these  scans  offer  insight  into  the 
actual  composition  of  the  sample,  as  well  as  give  quantitative  information  that  can  be  used  in  the  context  of  forward 
models.  The  approximation  method  shown  here  has  potential  to  expedite  these  forward  model  solutions  by  providing 
the  ECT  information  over  the  microstructure.  This  bodes  well  for  future  inversion  schemes  that  may  take  advantage 
of  the  physical  intuition  of  anisotropic  material  noise.  Incorporation  of  this  into  parameter  estimation  routines  remains 
as  future  work  First  and  foremost,  a  very  thorough  investigation  of  the  validity  of  the  assumptions  made  must  be 
undertaken  and  the  approximation  method  must  be  validated  with  real  experimental  data. 


Coil 

=  R\<j 

ct2  =  R2o 

FIGURE  9.  Image  of  the  geometry  analyzed  with  the  FEM  model  as  well  as  the  approximation  method. 
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FIGURE  10.  Results  from  the  FEM  model  is  represented  by  dots,  and  results  from  the  synthetic  aperture  routine  are  shown 

as  a  line. 
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